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Abstract 

We present a theoretical study of the spectral interferences in the spontaneous Raman scattering 
spectra of major combustion products in 30-atm fuel-rich hydrogen-air flames. An effective 
methodology is introduced to choose an appropriate line-shape model for simulating Raman 
spectra in high-pressure combustion environments. The Voigt profile with the additive 
approximation assumption was found to provide a reasonable model of the spectral line shape for 
the present analysis. The rotational/vibrational Raman spectra of H 2 , N 2 , and H 2 0 were 
calculated using an anharmonic-oscillator model using the latest collisional broadening 
coefficients. The calculated spectra were validated with data obtained in a 1 0-atm fuel-rich H 2 - 
air flame and showed excellent agreement. Our quantitative spectral analysis for equivalence 
ratios ranging from 1.5 to 5.0 revealed substantial amounts of spectral cross-talk between the 
rotational H 2 lines and the N 2 ONQ- branch; and between the vibrational H 2 0(0,3) line and the 
vibrational H 2 0 spectrum. We also address the temperature dependence of the spectral cross-talk 
and extend our analysis to include a cross-talk compensation technique that removes the 
interference arising from the H 2 Raman spectra onto the N 2 , or H 2 0 spectra. 
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1. Introduction 

Quantitative measurements of species concentration and temperature in elevated pressure 
combustion environments are critical to the validation of computer models used in the 
development of gas turbine and internal combustion engines. Laser-excited Spontaneous 
Raman scattering (SRS) spectroscopy is perhaps the only practical technique that 
provides spatially and temporally resolved multi-species information in a simultaneous or 
single-shot fashion with a high level of accuracy [1-3]. For high pressure combustion 
studies, SRS has the added benefit of providing stronger signal levels than those 
performed at atmospheric pressure, because the signal scales linearly with the number 
density, and hence pressure of the molecules being probed. Furthermore, SRS 
measurements are relatively simple, less expensive, and more versatile when compared to 
coherent anti-Stokes Raman spectroscopy (CARS), or stimulated Raman spectroscopy. 
These facts have motivated a recent renewed interest in fundamental studies that apply 
SRS diagnostics in high-pressure combustion [4-5]. 

Although the use of SRS diagnostics for quantitative studies in flame studies is well 
known, special attention to spectral interferences between Raman scattering from 
different molecules will still require further study. This requirement is substantiated by 
the fact that several earlier studies have indicated that spectral interferences can be a 
significant problem for Raman diagnostics in flames [2-3, 6-8]. In fact, the use of a 
detailed spectral analysis that includes the spectral interferences or ‘cross-talk’ between 
different vibrational/rotational Raman bands of different molecular species, is the basis 
for several common SRS techniques including thermometry via Stokes/anti-Stokes ratio 
[9], and integrated area measurements of vibrational Raman bands for species 
concentrations [10]. Without considering the spectral interferences among the different 
Raman bands for multiple species, the accuracy of temperature/species measurements is 
insufficient for combustion model validation. 

The effect of spectral interferences on the SRS signals has been effectively 
implemented using the well-known cross-talk matrix calibration procedure [2-3,6-8]. 

The off-diagonal terms of the calibration matrix generally represent the Raman scattering 
crosstalk of one species onto another. The acquisition and development of an accurate 
cross-talk calibration matrix is the key factor in obtaining quantitative species 
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concentration and temperature measurements in flames. Unfortunately, the Raman cross- 
talk calibration data obtained from these one-of-a-kind setups to date is not easily 
convertible into a transferable standard that is suitable for use by other researchers, 
particularly those who deal with high-pressure flames. In order to create a standardized 
cross-talk calibration matrix with absolute numbers, a comprehensive and accurate study 
of the spectral interference among Raman scattering in flames of various conditions (fuel- 
air composition, equivalence ratio, and temperature) over a range of pressure is still 
required. 

The quantitative treatment of the spectral interference in flames is not an easy task 
since the interference is a molecule-dependent function of temperature, species 
concentration, and chemical composition. For instance, rotational H 2 Raman lines at 
rotational quantum number, J=9toJ = 11 around the 2300 cm" 1 Raman shift in fuel- 
rich hydrogen flames may become strong enough to interfere with the N 2 ^-branch 
spectrum (near 2330 cm" 1 ). In particular, the rotational J— 10 Stokes line of H 2 is 
completely overwhelmed by the N 2 ^-branch. In order to facilitate the quantitative 
analysis of these interferences, a theoretical means of predicting the SRS spectra is 
essential. The quantitative analysis of the interferences are even further complicated in 
high-pressure flames due to possible pressure broadening or collisional line-narrowing 
effects [1 1]. 

In this paper, we investigate the spectral interferences between combustion 
molecules in high-pressure, fuel-rich, hydrogen-air flames using a spectral simulation of 
the SRS process. We focus our attention in this paper to hydrogen-air flames for the 
following reasons: chemically, a FU-air flame provides the simplest combination of 
combustion products in order to facilitate the spectral analysis, yet retains spectral 
features that are important to more complicated fiiel-oxidizer systems that will be studied 
at a later time; research data on a Fh-air flames are important in high-speed propulsion 
systems such as rocket combustors and pulse detonation engines. Fuel-rich hydrogen 
flames are studied here because they are typically used in rocket engines for reasons of 
flame stability and maximum heat-release. Furthermore, the interference between H 2 and 
other molecules becomes important in fuel-rich flames where H 2 concentration is high. 
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We first introduce an analytical methodology to validate an appropriate line-shape 
model for simulating Raman spectra in high-pressure flames and compare it with 
experimental data to verify our theoretical model for simulating Raman spectra. We then 
model the rotational/vibrational Raman spectra of H 2 , N 2 , and H 2 O using an anharmonic- 
oscillator model with a line-shape model that incorporates the latest collisional 
broadening coefficient data. Finally, we theoretically investigate the spectral behavior of 
Raman 0-, S-, and g-branches of the species in hypothetical high-pressure hydrogen 
flames over a range of equivalence ratios. 

2. Modeling of Raman Spectra 
2-1. Raman Cross-Section 

The theory of Raman scattering in diatomic and polyatomic molecules is well established 
[12]. According to Placzek’s polarization theory of the Raman effect, the corresponding 
differential Raman scattering cross-section per molecule at thermal equilibrium is given 
by [12] 


r dcr'' 

9 

M 

K dCl j 


v £ 0 ; 


(v 0 ±v*) g 


4 exp(-£, / kT) ^ 2 


Q(T) 




( 1 ) 


Here, the superscript Vindicates the relative angle of the planes of polarization between 
scattered and incident light; the subscripts of i and / represent the initial and final states 
respectively; is the vacuum permittivity; is the incident light (laser) frequency; + 
and - signs refer to the anti-Stokes and Stokes Raman scattering respectively; v k is the 
molecular frequency associated with the rotation-vibration, or the pure rotational 
transition from i to /; g is the total degeneracy; Et is the energy of the initial state; k is the 
Boltzmann constant; T is the gas temperature; Q(T) is the total partition function of the 
molecule at temperature T; p 2 ,0) is the molecule-dependent temperature- 
independent invariant that is a function of the mean value of the derived polarizability 
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tensor a 2 and the anisotropy tensor /T 2 . The values for d^o 2 , p 2 ,0) are tabulated by Long 
[12] for every branch of vibrational and pure rotational scattering for the molecules of 
interest; an dflj) is the correction factor arising from the rotation-vibration coupling [13]. 
Note that Eq. (1) should be multiplied by the factors (v + 1)Z> V 2 and vb v 2 for vibrational 
Stokes and anti-Stokes scattering of a diatomic molecule, respectively. Here, v is the 
vibrational quantum number and b v is defined as h/Sn c Vk where h is the Plank constant, 
and c is the speed of light. 

The Raman line frequency is given by v* = ( Ef - E t )/hc where Ef is the energy of the 
final state. The energy in a rovibrational state (v, J) is expressed as the sum of the 
vibrational energy term, G(v) and the rotational energy term, F(v, J) 

E(v,J) = hc[G(v) + F(v,J)]. (2) 


According to the anharmonic oscillator model, the vibrational energy term G(v ) of a 
diatomic molecule is given by [14] 


G(v) = co e 


( n 
l 2 ) 


v + — 


+ a e y e 




v + — 


(3) 


where o) e , co e x e , and co e y e are the molecular constants and are tabulated by Huber and 
Hertzberg [15]. According to the non-rigid rotator model that accounts for the influence 
of the centrifugal distortion, the rotational term for the *2 (singlet) state of diatomic 
molecules such as N 2 or H 2 is given by 

F(y,J) = B V J(J + 1) - D, J 1 (J + 1) 2 + H, J 3 (J + 1) 3 (4) 


where B, = B e - a e (v + 14) + y,(v + 14) 2 - rj e (v + !4) 3 , D V = D„ + fi e (v + 14) + ^ e (v + !4) 2 , 
and H v ~ H e - Q(v + Vi). The terms B e , D e , H e , a e , , y e , Tj e , and are the molecular- 
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dependent spectroscopic constants, which are tabulated for each molecule of interest by 
Huber and Hertzberg [1 5]. Thus for a Stokes branch (Av = v/- v, = +1), the Raman line 
frequency is given by 


v k( v t> J <) 


£(v, +l,J / )-£(v„J,) 
he 


( 5 ) 


The partition function Q(T) for a rovibronic state may be written as the sum of the 
rotational and the vibrational partition functions 


Q(T) = Q vlb (T) + Q rol (T) 


( 6 ) 


where 


Q vlh (T)= Xexpl 


v=0 


- hcGiv ) 
kT 


( 7 ) 


Q^T)=n s x2j^J^^) 

v=oy=o \ kl ) 


( 8 ) 


Here g s is the nuclear spin statistical weight and the term (2 / + 1) is the rotational 
degeneracy. 

For calculating a pure rotational Raman spectrum (e.g. rotational H 2 S-branch), Eqs. 
(2), (5), and (6) are replaced by E(J) = hcF(J) where F(J) - F( 0 , J), Vk(Ji) = - 

E(J,)]/hc, and Q(T) ~ Q ro t(T), respectively. 

The correction factor j{J) is given by 

f(J,) = 1 + 4 X -'(B, /« J 2 [j, (j, +\)+J f {J f + 1)] , (9) 
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where %~ l is the polarizability-related spectroscopic term. The values of % x for N 2 , and 
H 2 have been found to be 2.5 [13]. The factor j(J) is substantial for N 2 at temperature 
above 1000 K and is significant for H 2 because of its large rotational constant [13]. 

Polarizability Tensor Components 

According to Placzek-Teller theory, the space-averaged O (a 2 , fi 2 ,0) terms for the Raman 
Q-, O-, and ^-branches for scattered light that is polarized along the direction of the 
incident light polarization are given by the following equations, respectively [12,16-17]: 

®(«,/S) Q =| (v'.j l a lr. + i,j>f I Alrv.u>| 2 . 0°) 

S j-ij (¥rj P I Y v +\,j~i ) ’ (11) 


<D(a,/?) s — ^S j+ 2 j(y/ v j I P I Wv+xj+i) 


( 12 ) 


where \j/ v j is the vibrational wavefunction for a particular v-J state; a is the isotropic part 
of the polarizability, and /? is the anisotropic part of the polarizability. Sj,r are the 
Placzek-Teller coefficients for diatomic molecules given by [1,12] 


J(J + 1) 

(2 J - l)(2J + 3) 


(^-branch), 


(13) 


V.J-2 


3J(J - 1) 

2(2 J + 1)(2 J - 1) 
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Evaluation of Eqs. (10-12) based on the reasonable approximation for vibrational 
transitions with Av = 1 yields [16] 


<D(a,/?) Q =(B>J 


' a ' 2 +—Sj,p ' 2 

V 45 JJ 


(16) 


<^{a,P) 0 ={Bjco c )—S J 
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(17) 


0(a,p) s =(Bjo >t )^S J ^ J p' 2 . (18) 

Substituting Eqs. (16-18) into Eq. (1) and evaluating the degeneracy terms we obtain the 
following relations: 
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( 21 ) 


where the superscript ZZ indicates the scattered light polarized along the direction of the 
incident light polarization; q is a molecule-independent physical constant; g s is the 
nuclear spin statistical weight (g* = 1 for even J; g s = 3 for odd J ); (2J+1) is the rotational 
degeneracy. The terms a' and / 3 ' are the isotropic (mean) and anisotropic invariants of 
the static dipole polarizability derivatives, which are derived from the expansion in the 
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dimensionless vibrational coordinate £= (r-r e )/r e (where r is the intemuclear distance and 
r e is the equilibrium intemuclear distance) [16-17], 


a = a 0 +a'<^ + — a"^ 1 + ..., 


( 22 ) 




(23) 


with 


a' = da I84\ ( , 0 , 

(24) 

P' = dp!dp | { , 0 . 

(25) 


H2O Raman Model 

The Raman scattering spectrum of water vapor appears over the range of frequency from 
3400 to 4130 cm' 1 , and results mainly from the OH-stretch vibrational mode. Generally, 
the temperatures encountered in combustion environments makes the V\ band (100)-(000) 
transition dominant. For the H 2 O Raman spectrum, the classical anharmonic-oscillator 
model is not suitable for calculating the vibrational energy because of the asymmetric top 
nature of the water molecule that shows a relatively complicated structure. Fortunately, 
the vibrational energies in the vibrational ground state (000), and in the various excited 
states such as the (100), and (001) are known from the analysis of a high-resolution 
infrared spectrum and such database is tabulated in HITRAN [18]. The wavefunctions 
for those states were also calculated based on the rotational Hamiltonian by Flaud et al 
[19] that includes the Fermi resonance between the (100) and (020) states, and the 
Coriolis-like anharmonic coupling between these two states and the (001) state. Based on 
the wavefunctions and the polarizability derivatives, Avila et al [20] calculated the 
vibrational Raman cross-sections of H 2 0 in the OH-stretching region with great accuracy 
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up to about 1000 K. In the current report, we then used the results on the Raman line 
frequency and cross-sections which are tabulated in Ref. 20 to simulate the H 2 O Raman 
spectra in our high-pressure flames. 

The partition function of H 2 O was calculated using the well known polynomial 
expression [18,21] 

Q(T ) = a + bT + cT 2 + dT 3 , (26) 

where coefficients a , b, c, and d are tabulated over the range of temperature up to 3005 K 
by Gamache [21]. 

It should be noted that the resulting calculated H 2 O Raman spectrum does not 
include the hot bands of the (1 10)-(010) and (01 1)-(010) transitions that become 
important over 1500 K as well as the additional hot bands from the (020)-(000), (100)- 
(000), and (001)-(000), which would appear at higher temperature as Avila et al [20] give 
caution. These hot bands generally “stretch” a H 2 O spectrum out to the shorter 
wavelength, thus they form a longer blue-side tail. According to Avila et al [20], the 
intensities of these hot bands amount to over 1 0% of that of the main bands so that the 
present calculation for H 2 0, likely, also has errors of this level at higher temperature 
conditions. 


2-2. Line-Shape 

Pressure-Broadening Coefficient 

Collisional broadening coefficients of molecules are a dominant factor for the line-shape 
modeling. Recently collisional broadening coefficients for H 2 with perturber molecules 
have been proposed by advanced line-shape calculations and those numbers over a range 
of temperature and pressure have been successfully applied to CARS thermometry in 
high-pressure flames [22-24]. Bergmann et al [22] obtained J-dependent H 2 collisional 
broadening coefficients for higher temperature from 1800 to 2200 K for J- 0 to 7 in a 
fuel-rich CfL-air flame up to 40 bar. Hussong et al [23] provided an effective, J- 
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dependent polynomial equation for the H 2 collisional broadening coefficient including the 
inhomogeneous contribution (nonlinear dependence on the perturber concentration) in the 
H 2 -N 2 binary system: 

7j(T,C) = {ff)j +(.r*' t )jT + (a t ) J C + (a 2 ) J C 2 +(a i ) J CT + (a A ) J C 1 T, (27) 


where y] [cm'Vatm] is the total broadening coefficient as a function of temperature, T 

and concentration of perturber N 2 , C. Here C is defined by the N 2 fraction relative to H 2 
concentration under the assumption of the H 2 -N 2 binary system. All other parameters, 
(yf ), , (y !e,f )j , (a,)j (i = 1 to 4) are tabulated for J~ 0 to 5 in Ref 23. 

Since the influence of H 2 O as an additional collider plays an important role at higher 
temperature in terms of speed-inhomogeneous effects as Hussong et al [23-24] 
mentioned in their reports, a ternary system H 2 -N 2 -H 2 O should be considered for better 
accuracy over the range of flame temperature. Unfortunately, the collisional broadening 
coefficients of ternary mixture including H 2 O at high temperature have not yet been 
reported in the literature. We thus calculated the H 2 broadening coefficient for each 
flame condition using Eq. (27) in this study. Furthermore, in an effort to simplify the 
effects of temperature and 7-number dependence on the collisional broadening for 
hydrogen, we use the broadening coefficient of J= 3 for all lines in the H 2 spectrum. 
Note that the above-mentioned assumption regarding the binary system and the single- 
value coefficient has been previously found to be both appropriate and practical from a 
diagnostics standpoint [25]. 

There is an alternative expression for the collisional broadening coefficient, which 
has been proposed in the HITRAN database [18]. Molecular spectroscopic information 
(including the collisional broadening coefficients) of major combustion-related species 
such as N 2 or H 2 O are available in the database while there is no information regarding 
H 2 , due to the fact that it is not an infrared active molecule. We then used HITRAN [18] 
to find spectroscopic information for the collisional broadening coefficients of N 2 and 
H 2 O. In the database, the collisional broadened halfwidth is written as 
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~0 /nn -n HITRAN 

Yj (T,t) 


^ j [r cn . T 0 )C + r ,M (P»Jo )(1 - C)] 


(28) 


where f\PoJo\ and f cl \PoJ o) are the air-broadened, and the self-broadened width at 
the reference temperature, To - 296 K and the reference pressure, Po = 1 atm, and n is the 
coefficient of temperature dependence. As for N 2 , only f d \Po,To) was considered and it 
was assumed to a single value of 0.048 cm'Vatm. This value was derived from HITRAN 
by calculating an intensity-weighted, averaged halfwidth of all N 2 ^-branch lines from J 
= 1 to 40 in the vibrational 1-0 transition, which is the dominant transition at room 
temperature. For the H 2 O lines, we used the widths given by f x \Po,To) = 0.0095 and 
f e[f (Po,To) = 0.0923, which were also derived from the HITRAN database by calculating 
an average halfwidth of all lines in the vibrational transition (000)-(100) between 3400 to 
3710 cm' 1 at room temperature. A value of n = 0.5 was used for N 2 according to Ref. 18 
and n = 0.67 was found reasonable for H 2 O after calculating an averaged n value for the 
(000)-(100) transition. 


Process of Line-Shape Model Selection 

Spectral line shapes of Raman scattering spectra at atmospheric pressure are typically 
modeled by a Gaussian profile, and those at higher pressures are usually modeled by a 
Voigt profile which is the convolution of Gaussian and Lorentzian collisional broadening 
distributions. However, the Voigt profile is based on the impact approximation, so it may 
not be valid for modeling spectral shapes at very high pressures due to the break down of 
the approximations in the far wings of spectral lines [27], In some cases that require 
higher accuracy, such as in a laser absorption technique, motional ly-narro wed Galatry 
profiles [26] could be also employed. These line-shape models were developed for a 
single isolated transition. In practice however, multiple overlapping transitions are 
common and the effect of their overlap may become important at high pressures. Thus, 
consideration of the line-narrowing due to the interaction between different transitions, 
the so-called line-mixing effect, may be required to properly model the Raman spectra for 
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quantitative measurements. The line-mixing effect in nonlinear Raman processes of pure 
gases at high pressure has been previously studied [1 1,28] but not in flames. 

Based on the above points, we see the need to develop a methodology to choose an 
appropriate high-pressure line-shape model for a particular species, and to evaluate the 
limitations of these models. Nagali [29] discussed the applicability of the major line- 
shape models and introduced an effective approach for choosing an appropriate model for 
neutral species within the temperature and pressure range of interest. In this study, we 
follow the the approach described by Nagali [29] but include, for the first time, N 2 
concentration-dependent collisional broadening coefficients to the analysis. This permits 
the determination of an appropriate collisional broadening model for simulating Raman 
spectra in high-pressure flames. 

We consider three criteria for modeling the spectral lineshapes of Raman scattering 
spectra: 1) Motional narrowing (in which the Galatry profile is considered); 2) Impact 
approximation (in which the Voigt profile is valid); and 3) Line mixing (in which the 
additive approximation is not valid). To define the effectiveness and limitation of each 
criterion in terms of the range of temperature and pressure, we first define a collisional 
scale factor, i.e. the molecular collision duration, t c0 u [s]. When the interaction length is 
set to the radius of the collision-broadening cross-section, t co \\ can be expressed as 


T,C) 


1 

C' Av° 

V mol V 



(29) 


where Av° [cm' 1 ] is the collisional broadening fullwidth which is twofold of the total 

broadening coefficient at room temperature and atmospheric pressure, c [m/s] is the 

speed of light, TV, [cm' 3 ] is the number density of the molecule of species j; and v mo i [m/s] 
is the mean speed of the gas molecule which is approximately equal to another measure 
of the molecular speed, root-mean- square speed [30] 
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v mol (r) = 


3 RT 
V M 


(30) 


where ^ [J K 1 mol *] is the gas constant, M[g mol' 1 ] is the molecular mass. 

The wavelength region over which the impact approximation (collisional broadening 
due to intermolecular collisions which can be described by Lorentzian profile) is valid 
can be determined using the expression [31] 

2^„||Av imp «l , (31) 

where Av j mp is the spectral region in which the impact approximation is accurate. Thus, 
in the impact approximation region, zlvjmp can be written as 


&v imp (T,C) = 


0.1 


(32) 


Here we use a factor of 0.1, which is a smaller value than that suggested by Nagli (0.2) 
[29]. However, the exact value for this quantity is still open to interpretation for future 
studies. Following Nagali [29], if we require that at least 95 % of the area under the line 
shape be accurately described by the impact approximation, the criterion that the impact 
approximation is valid is then given by 

~2~ ~ 0-079Av imp . (33) 


Here, Av c [cm' 1 ] is the collisional broadening fullwidth given by 


Av c (T,C,P) = 2r;P , 


(34) 
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where P [atm] is the total gas pressure. In a regime above the boundary defined by Eq. 
(33), a line shape may still be modeled by Voigt function, but a breakdown at the far 
wings will likely be observed. 

The so-called line-mixing effects caused by the interaction between different 
transitions become important when the frequency separation of collisionally-coupled 
transitions, relative to consecutive rotational levels, is smaller than, or of the same 
magnitude as the collisional width [31]. That is, if neighboring, interacting transitions Vk\ 
and via are such that 


W 


*i 




(35) 


then line-mixing effects should be considered. Here, \v kX - v k2 | is the separation distance 

in frequency between the interacting transitions, and A v c is the averaged collisional 
fullwidth of the two different transitions of interest. Since we have to consider multiple 
lines in Raman spectra we used the averaged line separation, \v kl -v k2 \ = v^p between 

each consecutive rotational/rovibrational line for practical convenience. Thus, the 
boundary of the line-mixing can be written as 


v 


< 

sep 


2 


(36) 


Here, values for v sep were found to be 200.0, 19.8, 3.6, and 0.73 cm' 1 for the rot-U 2 (J=0 
to 10), vib- H 2 (J- 0 to 17 in vibrational 1-0 transition), H 2 O (000-100 transition), and N 2 
(g-branch, 7= 1 to 40 in vibrational 1-0 transition), respectively. 

Because of interrupted motions caused by the collisional scattering of molecules, a 
velocity-changing collision results in a line shape with a sub-Doppler linewidth, which 
involves line-narrowing effects [26]. If velocity-perturbing collisions are assumed to be 
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independent of state-perturbing collisions, then the Galatry profile (“soft” collision 
model) can be used to model the motional narrowing [29]. When the relationship given 
by 


Av c >Av 


(37) 


is true, here Avd [cm' 1 ] is the full-width-half-maximum (FWHM) Doppler width, 
collisional broadening effects dominate the Doppler broadening so that the refinement of 
the motional narrowing effect is not necessary [29]. This is the third criterion for the 
line-shape model. Here, the Doppler width is given by 


r\ 

A v D ( T ) = — ^2 ln(Z)kT/M , 
c 


(38) 


where vq is the transition center frequency. Note that the motional narrowing is different 
from the collisional narrowing in its overall shape due to the line-mixing effect at high 
pressure. 

From Eqs. (29-34 & 36-38), the following three boundaries are formed and they 
provide a way to choose the appropriate line-shape model for H 2 , N 2 , and H 2 O gases over 
a range of pressure, temperature, and species concentration: 


P = 


Ak„/2 , 
fj(C,T ) ' 


Galatry boundary - A 


(39) 


0.079 A v- 

P = — : Voigt boundary - B 


(40) 
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P = - Q sep — : Line-mixing boundary - C 

y j (c 9 


(41) 


With these boundaries, the region with Pf] <Av D /2 indicates the neccesity of 
considering the motional narrowing effect, i.e. use of the Galatry profile; that with 
Av d /2< Pf] < 0.079A v imp indicates the use of the Voigt profile for entire line (additive 

approximation is valid); that with 0.079A v jmp < Pf ] < v sep supports the use of the Voigt 
profile, with the caveat, that there is a possible breakdown in the far-wing (additive 
approximation is likely valid); and that with Pf] > v sep indicates the neccesity of 
considering the line mixing effect. 

2-3. Target Flames 

In order to analyze spectral interference in Raman scattering in high-pressure combustion 
environments, we generate thermodynamic properties for five different flames under the 
assumption of complete combustion and chemical equilibrium. Chemical equilibrium 
compositions of several different hydrogen-air flames, which hereafter we call “target 
flames”, were then calculated by the CEA code [32] for assigned equivalence ratios, 
pressure (30 atm), and an initial gas temperature (300 K). Table 1 shows the calculated 
mole fractions of H 2 , N 2 , and H 2 O and the adiabatic temperatures. The equivalence ratios 
were varied from 1 .5 to 5.0, with the resulting temperatures varying from 1402 K to 2262 
K. First, the calculated temperatures, pressure, and N 2 mole fractions as perturber 
concentrations were used as input parameters to calculate the collisional broadening 
coefficients and the line-shape diagrams. The calculated temperatures, pressure, and 
number densities of H 2 , N 2 , and H 2 O (derived from the mole fractions and temperatures) 
were then used for the Raman spectral simulations. 


3. Results and Discussion 
3-1. Choice of Line-Shape 
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To choose the best line-shape function to model Raman spectra in our target flames over 
a range of pressure, temperature, and perturber concentration we introduce a three- 
dimensional diagram that is consists of the boundaries defined in Section 2-2. Figure 1 
shows the diagram designed for the S(3) line of the pure rotational H 2 (hereafter indicated 
rot- H 2 ) spectrum ( Vo = 1035 cm' 1 ). This diagram directs one to choose an appropriate 
line-shape model for spontaneous Raman spectra simulation in any flame of interest with 
a prior knowledge of pressure, temperature, and N 2 concentration. The locations of the 
five target flames (labeled case # 1 to # 5 in Table 1) are marked in the diagram to 
determine which line-shape model is the most appropriate. All five cases here are 
between boundary A and B as shown in Fig. 1, which indicates the entire line shape can 
be modeled by the Voigt profile with the additive approximation. Note that the same is 
true for other major rot-Wi lines of higher J number while the Galatry boundary 
(boundary A) slightly moves toward higher pressure depending on the transition center 
frequency Vo. 

Figure 2 shows a three-dimensional diagram for the rovibrational H 2 (hereafter 
indicated as vib- H 2 ) spectrum with points representing the target flames. Here the 
fundamental frequency of the W 6 -H 2 spectrum (4160 cm’ 1 ) was assigned for the value of 
vo in the calculation of the diagram. In this figure, the line-mixing boundary (boundary 
C) appears because the averaged line separation (v sep ) of the vib- H 2 spectrum is 19.8 
cm' 1 , which is much smaller than rotational separation of the rot-. H 2 spectrum by 180 
cm' 1 . All the target flames in Fig. 2 are also between boundary A and B, the Voigt 
profile regime. 

Figure 3a and 3b show the line-shape diagrams for the vibrational (2-branches of N 2 
and H 2 O spectra with points representing the target flames. The fundamental vibrational 
Raman frequencies of N 2 and H 2 0 (2331 and 3657 cm' 1 , respectively) were assigned to 
Vo in the calculations of the diagrams. Here, only the self-broadening is considered for 
N 2 , and the N 2 fraction relative to H 2 0 concentration is constant, i.e. C = 0.65. 
Accordingly, the diagrams have only two dimensions without a dimension of perturber 
N 2 concentration. As shown in Fig. 3a, all the target flames are located between 
boundary B and C, so the line shape for the N 2 (7-branch can be appropriately modeled 
by the Voigt function with the assumption of the additive approximation in the near-wing 
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region, but their line shape could break down in the far-wing regions. The impact and 
additive approximation are, for the most part, appropriate for the entire region of the H 2 O 
spectrum in the target flames according to Fig. 3b. 

In general, the Voigt profile was found to be a reasonable model for simulating the 
Raman spectrum over the range of pressure up to 30 atm, and for temperatures from 
approximately 1500 K to 2500 K. While there may be a slight effect of breakdown at the 
far wings of the N 2 ^-branch, the effect is small and is ignored in favor of the simplicity 
of using the Voigt function for the entire region of interest. This consequence is very 
useful for reducing the computational requirements of the spectral simulation and 
analysis. Hence, we used the Voigt profile for modeling the line shape of Raman spectra 
in the target flames. In the sections below, we will empirically validate the use of the 
Voigt profile for simulating Raman scattering spectra via comparisons between 
calculated spectra and experimental data. 

In calculations of the Voigt-function a and x parameter [33] we used the overall 
convolution Gaussian width, Avq [cm' 1 ] given by 

A v G (v, J) = ^Av d 2 +Av/ spec 2 . (42) 


Here A v spec is the instrumental function, which is the spectral resolution of a 
spectrograph. Finally, the profile of the Raman spectrum at temperature (7), and pressure 
( P ) was calculated by the integration (via summation) of all Voigt-shaped lines based on 
the additive approximation assumption. 

Regarding the temperature dependence of the collisional broadening coefficient there 
are at least two alternative expressions (CARS-, or HITRAN-expression) as we discussed 
above. We are not aware of a published report that indicates whether or not these two 
expressions make any difference to the line-shape diagram. In Fig. 4, we compare the 
line-shape diagrams derived from the two different expressions for the broadening 
coefficient. As shown in Fig. 4, the boundary A, and B based on the two expressions, 
which have been validated by the different spectroscopic experiments and calculations, 
agree very well when the HITRAN-based temperature dependence coefficient, n was set 
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to 0.15. This agreement indicates the both expressions are comparable, and that they can 
be applied to spontaneous Raman scattering applications, even though we used a CARS- 
derived expression for H 2 , and a HITRAN-derived expression for N 2 , and H 2 O in this 
paper. 

3-2. Validation of the Raman Model 

We validate our theoretical models for the Raman scattering intensity and line shapes in 
high-pressure environments by comparison with experimental data measured in high- 
pressure flames. We used a specially designed high-pressure gaseous burner that 
provides a precise and stable source of hot combustion products in a high-pressure 
environment. The burner provides a homogeneous stream of combustion products via a 
tightly-spaced array of multiple fuel-rich H 2 -air premixed flamelets (10 atm). In this 
particular premixed burner design, the 7 x 7 array of premixed flame tubes (1.1 mm dia.) 
had a center-to-center spacing of 2.60 mm. Spontaneous Raman scattering from H 2 , N 2 , 
and H 2 O were then measured at a reduced-temperature, post-flame zone using our SRS 
diagnostic system which consists of a pulsed Nd:YAG laser, a laser pulse-stretcher, a 
high-speed electro-mechanical shutter, and a spectrometer fitted with a non-image- 
intensified CCD camera. For the simulations, we used an instrumental resolution of v^ec 
= 35 cm' 1 which corresponds to the spectral resolution of the spectrometer used in the 
measurement. It should be noted that our previous measurements with this high-pressure 
burner indicated that a condition close to the chemical equilibrium was reached over the 
range of equivalence ratio at the measurement location, however, the measured 
temperatures were below the adiabatic equilibrium temperatures by more than 100 K 
[34]. Further details on the burner and diagnostics have been described in our previous 
reports [25, 34-35]. 

For a quantitative simulation of a Raman spectrum, we empirically determine the 
factor («'//T) 2 , the squared ratio of the isotropic term to the anisotropic term of the 
polarizability derivative, and we verify the relative intensity of the ^-branch to the S/O- 
branches in vibrational N 2 or H 2 Raman spectrum. Note that calculated S'-/ 0-branches of 
N 2 , or H 2 Raman spectrum using theoretically determined values of a ' and /?' [16-17] do 
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not agree with our experimental data, thus we determined these values empirically as 
follows. We measured the factor {a' Ip') 2 for N 2 , or H 2 vibrational Raman spectrum in 
ambient-temperature air and compare it with a calculated spectrum at the same condition. 
As shown in Fig. 5, we found the best agreement in terms of the relative intensity of the 
(7-branch to the O-, and S-branch between data for the room air and a calculated N 2 
Raman spectrum at 300 K when using the factor for (a'/fi') set at a value of 0.25. The 
suitability of this value for N 2 was also substantiated by the good agreement between 
measured data and calculated values for a N 2 Raman spectrum in a 1 0-atm EE-air flames 
(<f> - 2.06) as shown in Fig. 6. We found that the factor {a' Ip') 2 for H 2 , when set at 0.35, 
provided the best agreement between a calculated O-branch intensity of the vib- H 2 
spectrum and a measured spectrum from a 1 0-atm H 2 -air flame {(/>- 3.0) as shown in Fig. 
7. 

We then compare the experimental Raman spectra measured in 1 0-atm fuel-rich EE- 
air flames with theoretical Raman spectra calculated at the temperatures determined 
experimentally by the line-intensity distribution of the rot- H 2 Raman spectrum [13,25]. 
Figure 6 shows the measured (^=2.06) and calculated (T = 1425 K) results of the rot- H 2 
and the vibrational N 2 Raman spectra. The transition frequencies (wavelength) and line 
shapes of calculated spectra show excellent agreement with experimental data for both 
species. As indicated in this Fig. 6, the N 2 O-branch is a possible spectral interference 
(cross-talk) with the S( 9) and S(8) lines of the rot-Ui band. Also, the rot- H 2 £(10) line, 
which is barely seen in this figure, is completely overwhelmed by the N 2 (7-branch. We 
will examine these spectral interferences in detail in the section below. 

Figure 7 shows the measured (^= 3.0) and calculated (T = 1 133 K) results of the 
vibrational H 2 O and the vib- H 2 Raman spectra. The calculated spectra for both species 
show excellent agreement with the experimental data in terms of the transition frequency, 
intensity and line shape. Note that we can clearly observe the spectral interference 
between the H 2 O spectrum and the vib- H 2 0(0,3) line around 657 nm. This interference 
is important because it affects the determination of the H 2 O number density in fuel-rich 
high-pressure flames, thus we will examine it in detail in the following section. 

Based on the comparisons between the spectral simulation and the experimental data, 
we have shown that the current theoretical models for Raman scattering line position, 
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intensity, and shape, are of sufficient accuracy to examine the effects of Raman spectral 
interferences amongst the major species of interest in high-pressure Th-air flames. 

Quantitative Comparison between Different Species 

In order to quantitatively analyze the contribution of spectral interferences in the overall 
Raman spectra for all species of interest, the relative scattering intensities of each 
different species has to be realized and properly implemented in the calculations. To do 
this, we define an empirical scaling factor for the scattering intensity of a rotational or 
vibrational Raman spectrum of species j relative to the N 2 Raman scattering intensity as 
follows: 




r n=a 




(43) 


where ms l^ p is the pixel intensity (CCD count per pixel) of a measured rotational or 

vibrational Raman spectrum of species j at equivalent ratio of (j> and pressure of P; n is 
the pixel number and a, b, c, or d indicate a corresponding pixel range for the summation 
of the spectrum of species j; Nj {(j), P,T ms ) is the calculated chemical-equilibrium number 

density of species j for (j), P, and the measured temperature T ms via the rot-Wi spectrum. 
The subscript, *N2’ indicates the molecular nitrogen species used as the reference. 

Based on the scaling factor Tj and the intensity-normalized calculated Raman 
spectrum, the scale-corrected intensity of a theoretical Raman spectrum for species j at 
pressure of P and temperature of T can be obtained by 


cal r j (P,T) = r j N J (</>,P) 


cd I J (P,T) 

Y^JJpjyl’ 


(44) 
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where N ;(</>, P) is the adiabatic number density of species j for ^and P; and cal Ij(P,T ) 

is the intensity of the calculated Raman spectrum at P and T. 

Based on Eq. (43) and the measured Raman scattering data obtained in 10-atm H 2 - 
air flames (Fig. 6, and Fig. 7), we were then able to determine the scaling factor for each 
species. Note that the instrumental function did not affect the resulting scale -corrected 
spectra, since the spectral intensity response of the entire optical system used for the 
measurements are calibrated in both intensity and wavelength using standard reference 
lamps. We then make the assumption that the combustion products have reached 
chemical equilibrium, and also with use of the scaling factors described above, and by 
using the parameters found in Table 1, we calculate the scale-corrected rotational or 
vibrational Raman spectra of N 2 , H 2 , and H 2 O in 30-atm fuel-rich FE-air flames. This 
finally enables a quantitative comparison of the spectral interferences between the 
different species of interest. 

3-3. Spectral Interference 

According to the experimental data shown in Figs. 6 and 7, there are two kinds of 
noticeable spectral interferences in Raman scattering in fuel-rich hydrogen flames: 
interference between the higher ./-number (8 <J< 11) rot- H 2 lines and the N 2 spectrum; 
and between the vib- H 2 0(0,3) line and the vibrational H 2 O spectrum. The behavior of 
Raman spectral interference is typically nonlinear since the rate of spectral interference is 
a function of both temperature and number density of at least two species in a spectral 
range of interest. For instance, in fuel -rich flames, rotational H 2 line intensities of 
transitions with higher J states are determined by a trade-off between the Boltzmann 
rotational distribution and the number density which is governed by the local 
stoichiometry. Thus, one of our questions is in which case does the Raman interference 
become more important: FE-rich, lower-temperature flames, or nearly stoichiometric, 
higher-temperature flames? In order to examine this question in a quantitative fashion, 
we simulated rotational and/or vibrational Raman spectra of H 2 , N 2 , and H 2 O in high- 
pressure hydrogen-air flames at different equivalence ratios using the theoretical 
formalism presented above. 
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Interference between rot-H2 and vib-N2 

Figure 8 shows the calculated Raman spectra of the rot- H 2 5-branch and the 
vibrational N 2 spectra in 30-atm H 2 -air flames from </>= 1.5 to 5.0. Occurrences of 
spectral interference among these spectra can be clearly seen in this series of spectral 
simulations. Our simulation indicates that it is particularly important to quantitatively 
analyze the interference that the rot- H 2 5(10) line gives to the N 2 (9-branch since the 
5(10) line is completely overlapped by the N 2 line. Figure 8 also confirms that there is an 
interference between the N 2 O-branch (the blue-side tail) and the lines from the rot- H 2 
5(8), or 5(9) lines; and that there is an interference between the N 2 5-branch (the red-side 
tail) and the rot-U 2 5(1 1) line. Note that the higher ./-number (> 9) line intensities of the 
rot- H 2 do not necessarily become stronger with increase of the equivalence ratio, while 
the lower J-number (< 7) line intensities apparently increase with the equivalence ratio. 
The interference between the ro/-H 2 and the N 2 spectrum indicates that further attention 
should be paid to the spectral shape of the N 2 (9-branch in hydrogen-rich environments. 
This effect is further exacerbated by the fact that the N 2 (9-branch gradually shrinks as 
temperature lowers due to the less population of higher vibrational states. 

In order to quantify and clarify the flame conditions in which the interference is 
significant, we define the cross-talk rate as the overlap rate between one spectrum and 
another. The cross-talk rate, K x is the ratio of the integrated spectral area of one spectrum 
(always smaller one) to that of another spectrum, which interferes with the former 
spectrum, for a certain range of wavelength. The spectral range of the integration for 
each case is shown in the figure captions. Figure 9 shows the cross-talk rate between the 
N 2 0-branch and the rot- H 2 5(8) line (hereafter indicated as /^ r N ‘ Hs8 ); that between the N 2 
0-branch and the rot- H 2 5(9) line (indicated as K r n ' Hs9 ); and that between the rot- H 2 5(10) 
line and the N 2 ^-branch (indicated as £ r Hsl0 " N ) over the range of equivalence ratios used 
in this study. A high amount of the cross-talk, A^ r N " Hs8 = 0.48 (48% overlap within the 
spectral range) is observed at (j) - 1 .5. £ r N ' Hs8 becomes smaller as (/) increases, yet x r N Hs8 
is more than 0.08 for (f) over 4.0. Similarly, the cross-talk Xr N ‘ Hs9 reaches more than 0.65 
(65% overlap) at (f)= 1.5. The higher amount of cross-talk arising from the “hot bands” 
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of the N 2 (9-branch also contributes to the interference. AT r N " Hs9 has a minimum value of 
0.3 (30% overlap) at (f) = 3.0 and stays over this value over the range of $ examined here. 

It is evident from Fig. 9 that the S(8) and S{ 9) lines should not be included in the 
temperature determination process via the rot-Wi distribution in fuel-rich flames because 
of the higher level of interference. 

In Fig. 9, the cross-talk rate ^ r N ' Hsl ° shows a different behavior than the previous 
cases for the S( 8) and S(9 ) lines: it varies from 0.023 (2.3% overlap) to its maximum 
value around 0.036 (3.6% overlap) at <f)~ 2.5. A cross-talk level of several percent is not 
insignificant, and may affect the accurate measurement of N 2 number density and/or N 2 
vibrational temperature, if this effect is not accounted for in the subsequent analysis. 

As an aid in the interpretation of Raman spectra with such an interference, we 
propose a practical method to compensate the cross-talk effect of ^ r N Hsl ° using a 
theoretical calculation with a priori knowledge of the temperature and an observed rot- H 2 
line intensity data. Since the intensity distribution of the rot- H 2 Raman transitions is 
governed only by temperature under assumption of the thermal equilibrium, the relative 
intensity ratio of any pair of lines in the rot-U 2 Raman spectrum can be theoretically 
derived based on a known temperature. Hence one can estimate the line intensity of the 
rot- H 2 S{ 1 0) from a measured intensity of, for instance, the rot- H 2 S( 6) line, which is 
generally strong enough to be measured from background noise and free from any major 
spectral inference. The theoretical Raman intensity ratio of the rot- H 2 5(10) line to the 
S(6 ) line, R r H sl0 ' s6 , was calculated with Eq. (21) over a range of temperature. The 
polynominal coefficients of the curve fitted to the theoretical plots of 7? r H Sl ° s6 are shown 
in Table 2. The thermodynamic temperature can be measured by the intensity 
distribution of the rot-U 2 Raman spectrum [13] within a few percent accuracy [25] and 
thus, a corresponding intensity ratio between the two lines at the temperature can be 
obtained via the polynominal function in Table 2. The “hidden” contribution of the rot- 
112 S(10) to the N 2 (9-branch would be then revealed by multiplying the obtained intensity 
ratio to the measured the S(6) line intensity in a Raman spectrum of interest. Finally, the 
obtained S( 10) line intensity would then be subtracted from a integrated intensity of the 
N 2 spectrum to find a corrected N 2 number density. 
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Note that the rot-H 2 S( 7) shown in Fig. 8 is not suitable for use in this cross-talk 
compensation process because it coincides with the common atomic Na emission line 
( 3s-3p transition) near 589 nm, which can be quite intense in industrial flames. For the 
other important cross-talk contributors to the vibrational N2 spectrum resulting from the 
S(9), or S(\ 1) lines in H2, a cross-talk calibration matrix formalism as previously 
mentioned in the Introduction section would have to be implemented in order to 
compensate for these effects. 


Interference between vib-H2 and vib-IfO 

Figure 10 shows the calculated vibrational Raman spectra of H2 and H2O in 30-atm 
H2-air flames over the range of equivalence ratio ranging from 1 .5 to 5.0. As the 
equivalence ratio increases (thus, the number density of hydrogen molecule increases), 
the intensity of the vib-Yii spectrum becomes stronger while the intensity of the H2O 
spectrum do not change much. As clearly seen in Fig. 11, the vib-Yh 0(0,3) line 
significantly interferes with H2O, particularly in highly Fb-rich flames. There is also 
interference around 661 nm by higher rovibrational transitions of the vib- H2 ^-branch. In 
order to clarify these interferences, we obtain the cross-talk rate between the vib- H2 and 
the H2O spectra (^ r vH W ) following the above-mentioned manner. The result over the 
range of equivalence ratios is shown in Fig. 11. The curve in Fig. 1 1 indicates that the 
cross-talk rate ^ r vH ‘ w has a monotonic and positive relation with equivalence ratio. 

Notice that the cross-talk rate x r vH ' w has a value of 0.08 (8% overlap) at (f) = 1.5 and 
reaches 0.18 (18% overlap) at <j)= 5.0. This level of cross-talk causes major errors in the 
measurement of H2O number density if not properly accounted for. 

Since the V/6-H2 0(0,3) line are mostly overwhelmed by the stronger H2O spectrum, 
it is generally hard to realize this interference. We then propose the same analytical 
compensation technique previously used for the cross-talk rate ^ r N Hsl ° to remove the 
contribution by the vib- H2 0(0,3) line to the H 2 0 spectrum. Using the a priori knowledge 
of temperature determined by the rot-Wi Raman transitions and a measured vib- H2 line 
intensity, a corresponding vib- H2 0(0,3) line intensity can be obtained and subtracted 
from the H2O spectral intensity data. Here, we suggest use of the V/6-H2 0(0,5) line for 
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the compensation process because it is the strongest line amongst the vib- H 2 O-branches 
and is relatively free from major spectral inferences. The theoretical Raman intensity 
ratio of the vib-Yij 0(0,3) line to the 0(0,5) line, i? V H° 3 ’ 05 , was calculated over a range of 
temperature. The polynominal coefficient of the fitting curve to the theoretical plots of 
i?vH° 3 ° 5 is shown in Table 2. The value calculated from the polynominal function in 
Table 2 with a known temperature and a vib- H 2 0(0,5) intensity can serve to compensate 
the H 2 O spectral intensity in fuel-rich flames. 

Another source of interference that arises from the blue-side tail of the vib- H 2 Q- 
branch around 661 nm, can cause a minor overestimation in measurements of H 2 O 
number density by less than -1.5% (e.g. at ^ — 5.0) after compensating for the vib-Vi .2 
0(0,3) contribution. There is not a simple, straightforward analytical expression for this 
interference because this portion of the vib-Hi g-branch contains multiple lines that 
originate from different rovibrational transitions, and its dependence on the spectral band 
factor (spectral resolution and/or pressure broadening) is complicated. Thus, the cross- 
talk matrix calibration method needs to be applied to take this interference into account. 


Temperature dependence of spectral interferences 

It is well known that the off-diagonal terms of the cross-talk calibration matrix 
represent spectral interference of one species onto another and may be expressed as a 
function of temperature [2-3, 6-7]. Although the temperature dependence of the cross-talk 
matrix has been typically assumed to be a polynomial relation [6-7], we are not aware of 
any studies that can substantiate this common assumption for high-pressure combustion 
environments. Figure 12 shows the cross-talk rate normalized by the number density of 
hydrogen (KjN m ) as a function of temperature regarding the N 2 ^-branch and H 2 0 

Raman spectra. By normalizing by the number density of H 2 , which is the only species 
whose Raman spectrum (either rotational or vibrational) interferes with the N 2 or H 2 O 
Raman spectrum in H 2 -Air flames, a pure dependence of the cross-talk on temperature 
can be extracted. From Fig. 12 we confirm that the temperature dependence of the H 2 - 
normalized cross-talk for nitrogen and water is well approximated by a monotonic, 
second-order (or less) polynomial function for the hydrogen flames in this study. 
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Note that the spectral database for the H 2 0 Raman scattering used here does not 
include a treatment of the hot bands, which is significant on the blue-side of the H 2 0 
spectrum above a temperature of 2000 K. As a result, the database tends to 
underestimate the spectral intensity by no less than 10% in high temperature regions. 
Thus, we have to be aware that an experimentally observed H 2 0 Raman spectrum will be 
more spread out towards the blue-side than predicted by our calculation. This hot-band 
behavior of H 2 0 could interfere with the vib- H 2 0(0,4) line in the equivalence ratio range 
from about 1.0 < ^ < 2.0, in high-temperature flames. According to our experimental 
data (with the same spectral resolution of the present simulation) in 10-atm H 2 -air flames 
over the wide range of equivalence ratio (0= 0.2 - 5.0) [34], the blue-side tail end of H 2 0 
spectrum touches the vib- H 2 0(0,4) line but barely interferes even at (f>= 1.5. Hence it is 
not necessary to include the vib-Wi 0(0,4) in the cross-talk calibration process for H 2 0 
Raman spectrum. 

4. Conclusions 

In order to improve and make spontaneous Raman scattering diagnostics more 
quantitative and accurate for the measurement of temperature and species in high- 
pressure combustion, we theoretically investigated the spectral interference of Raman 
spectra in fuel-rich, H 2 -air flames at pressures up to 30 atm. The rotational/vibrational 
Raman scattering of H 2 , N 2 , and H 2 0 were modeled with particular attention paid to the 
0-/iS-branches as well as the common ^-branch using an anharmonic oscillator model 
and the most recent spectral database available in literature. We applied an analytical 
methodology, which is based on a collision-scale analysis and molecular vibration, to 
validate an appropriate line-shape model for simulating the Raman spectrum in high- 
pressure flames. The line-shape diagrams for the species showed that the use of the 
Voigt profile with the additive approximation was appropriate and practical for the 
Raman spectral analysis at the target flame conditions used in this study. 

Our model was validated by direct comparisons between the calculated Raman 
spectra and measured data in 10-atm H 2 -air fuel-rich flames, which showed excellent 
agreement in every detail, including the 0-, and 5-branch intensities. The spectral 
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interferences amongst the Raman spectra of the species of interest were then calculated 
for the target high-pressure Fh-air flames, and were analyzed over the range an 
equivalence ratio range of <f> - 1.5 to 5.0, and temperature. Below is a summary of the 
findings of the quantitative analysis on the Raman spectral interference in 30-atm fuel- 
rich H 2 -air flames: 

(i) The spectral overlap (cross-talk rate) between the pure rotational H 2 £( 8) line and 
the N 2 (2-branch was as high as 48% at (f) = 1 .5. 

(ii) The spectral overlap between the pure rotational H 2 S( 9) line and the N 2 O-branch 
was 65% at (j> = 1.5 and about 30% at ^ = 3.0. 

(iii) The spectral overlap between the pure rotational H 2 *5(10) line and the N 2 (2-branch 
was 3.6% at ^ — 2.5. 

(iv) The spectral overlap between the vibrational H 2 0(0,3) line and the H 2 O spectrum 
was 18% at 5.0. 

In addition, we introduce an analytical technique using the Raman intensity ratio for 
a neighboring pair of the H 2 lines, to compensate for the cross-talk between the rot-Y^ 
£(10) line and the N 2 (2-branch; and that between the vib-Wj 0(0,3) line and the H 2 O 
spectrum. We also provide numerical coefficients for implementing this compensation 
technique. Finally, we demonstrated that the Raman cross-talk between rot- H 2 and N 2 , or 
vib- H 2 and H 2 O has a second-order (or less) polynomial relation with flame temperature; 
this information can then be used as an aid in selecting an appropriate functional form for 
calibrations in future studies. 

Since the collisional broadening effects on spectral line-shape are relatively minor at 
the modest spectral resolution considered here, the findings and techniques discussed in 
this report are applicable to both atmospheric and high-pressure (at least up to 30 bar) 
flames. 
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Table 1 Adiabatic temperatures and species compositions calculated for 
hydrogen-air combustion system using the CEA chemical equilibrium 
code using assigned pressures and equivalence ratios. 


Case # 

1 

2 

3 

4 

5 

P [atm] 

30 

30 

30 

30 

30 

<f> 

1.5 

2.0 

3.0 

4.0 

5.0 

T[ K] 

2262 

2065 

1770 

1560 

1402 

H 2 [Mole frac.] 

0.147 

0.257 

0.409 

0.509 

0.58 

H 2 O [Mole frac.] 

0.295 

0.257 

0.205 

0.17 

0.145 

N 2 [Mole frac.] 

0.55 

0.479 

0.381 

0.317 

0.271 
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Table 2 Polynomial coefficients of curves fitted to the theoretically calculated intensity 
ratios of lines for the rotational hydrogen S- branch and vibrational hydrogen Obranch 
as a function of temperature (K). The polynomial coefficients have the form of 
R(T) = aT A + bT 3 +cT 2 +dT + e. 


Coeff. 

a 

b 

c 

d 

e 

n sl0-s6 


-1.813-10' 11 

1.362-10’ 7 

-1.676-1 O' 4 

5.818-10' 2 

n o3-o5 
'HH 

4.089 10' 13 

-3.354-10' 9 

1.047T0 -5 

-1.508*10 2 

9.504 
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Fig. 1 Three-dimensional diagram used for choosing the proper line-shape 
model in high-pressure combustion environments for the S( 3) line of pure 
rotational H 2 Raman spectrum. Pressure is in log scale. White circles 
indicate the locations of the five target flames in shown in Table 1. 
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Fig. 2 Three-dimensional diagram used for choosing a proper line- 
shape model in high-pressure combustion environments for vibrational 
H 2 Raman spectrum. Pressure is in log scale. White circles indicate the 
locations of the five target flames shown in Table 1. 
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Fig. 3 Two-dimensional diagrams used for choosing a proper line- 
shape model in high-pressure combustion environments for (a) 
vibrational N 2 Raman ^-branch (self-broadening), and (b) vibrational 
H 2 0 Raman (C = 0.65). White circles indicate the locations of the five 
target flames shown in Table 1. 
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Fig. 4 Comparison of the boundaries in the diagram for the rotational 
H 2 S( 3) line with different pressure broadening coefficients based on the 
CARS-expression and the HITRAN-expression. Solid line : CARS- 
expression; dashed line : HITRAN-expression with temperature 
coefficient n - 0.15. 
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Fig. 5 Measured and calculated Raman spectra of vibrational N 2 in room- 
temperature air (300 K) at 10 atm used to determine the ratio of 
mean/anisotropic polarizability derivatives, a' //?' of nitrogen. Experiment 
(points): 200-shot averaged, 0.1 mm spectrometer slit. Calculation (line): v spec 
= 35 cm' 1 . 
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Fig. 6 Measured and calculated vibrational N 2 and rotational H 2 Raman 
spectra in a 10-atm H 2 -air flame. The spectra were calculated using the 
experimentally determined temperature via the rot- H 2 transitions, and was 
found to be 1425 K (v spec = 35 cm’ 1 ). The experimental data is represented 
by the points (^= 2.06, 200-shot averaged, 0.1 mm spectrometer slit). Solid 
line : rotational H 2 . Dashed line : vibrational N 2 . Two calculated lines are 
shown separately but not superimposed. 
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Fig. 7 Measured and calculated Raman spectra of H 2 0, and vibrational H 2 in 
a 10-atm H 2 -air flame. The spectra were calculated at an experimentally 
determined temperature obtained via the rot-H 2 transitions (not shown in this 
figure), and was found to be 1133 K ( v spec = 35 cm" 1 ). The experimental data 
is represented by points {(/>= 3.0, 200-shot averaged, 0.1 mm spectrometer 
slit). Solid line : vibrational H 2 . Dashed line", vibrational H 2 0. Two 
theoretical lines are shown separately but not superimposed. 
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Fig. 8 Calculated vibrational Raman spectra of H 2 0 and H 2 in 30-atm H 2 -air 
combustion at different equivalence ratios for the five flame cases shown in 
Table 1 with v spec - 35 cm-1. Dotted line is rot-H 2 , Dashed line is N 2 , and 
Solid line are superimposed (summation of the two spectra). Intensity is in 
arbitrary units but in the same scale for all equivalence ratios. 
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Fig. 9 Cross-talk rate (K r ) between N 2 O-branch and rot- H 2 lines in 30-atm 
H 2 -air combustion at different equivalence ratios, calculated from the spectra 
shown in Fig. 8 . Triangles are £ r N ' Hs8 ; Squares are K r N ' Hs9 ; and Circles are 
A^n-hsio. The spectral ranges for the determination of K t for S( 8), 5(9), and 
5(10) are 592.7 ~ 596.3, 598.1 ~ 602.0, 602.3 ~ 609.2 ran, respectively. 
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Fig. 10 Calculated Raman spectra of vibrational H 2 and H 2 0 in 30-atm H 2 -air 
combustion at different equivalence ratios for the five target flames in Table 1 with v spec 
= 35 cm* 1 . Dotted line is vib- H 2 , Dashed line is H 2 0, and Solid line is superimposed line 
(summation of the two spectra). Intensity is in arbitrary units but in the same scale for 
all equivalence ratios. Maximum values of vib- H 2 for the truncated peaks are included. 
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Fig. 1 1 Cross-talk rate (K r ) between vib-H 2 and H 2 0 spectrum in 30-atm H 2 - 
air combustion at different equivalence ratios, calculated from spectra shown 
in Fig. 10. The spectral range for the determination of K x is from 650.0 to 
664.0 nm. 
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Fig. 12 Temperature dependence of the cross-talk rate (K r ) in 30-atm H 2 -air 
combustion. Here the cross-talk rate is normalized by the adiabatic number 
density of H 2 for each condition to highlight the true temperature 
dependence. Circles are A^ r N ' Hsl °; Triangles are A" r vH ' w . Lines are second- 
order polynomial fit to the points. 


Kojima and Nguyen, submitted to JQSRT 




